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We study optimally doped Bi-2212 (Tc = 96 K) using femtosecond time- and angle-resolved 
photoelectron spectroscopy. Energy-resolved population lifetimes are extracted and compared with 
single-particle lifetimes measured by equilibrium photoemission. The population lifetimes deviate 
from the single-particle lifetimes in the low excitation limit by one to two orders of magnitude. 
Fundamental considerations of electron scattering unveil that these two lifetimes are in general 
distinct, yet for systems with only electron-phonon scattering they should converge in the low- 
temperature, low-fiuence limit. The qualitative disparity in our data, even in this limit, suggests 
that scattering channels beyond electron-phonon interactions play a significant role in the electron 
dynamics of cuprate superconductors. 

PACS numbers: 74.72.-h, 78.47.J-, 71.38.-k 


Electron lifetime is a central quantity in condensed 
matter theories [if. It determines macroscopic proper¬ 
ties such as electrical and thermal conductivities, and en¬ 
codes microscopic scattering mechanisms [H-Q. Reveal¬ 
ing dominant scattering channels in copper-oxide high- 
temperature superconductors (cuprates) will be key to 
understanding the complex interplay of orders underly¬ 
ing their phase diagram. 

Electron lifetimes in cuprates have been studied in 
both the energy and time domains. In the energy domain, 
angle-resolved photoelectron spectroscopy (ARPES) [3- 

and optical spectroscopy [13, El access the imagi¬ 
nary part of electron self energy ImS(e), which is con¬ 
nected to the single-particle lifetime 'rs(e) via ImE(e) = 
7i/(2Ts(e)). This lifetime describes the relaxation pro¬ 
cess of an excited single particle with energy e. On 
the other hand, time-resolved reflectivity (trR) [IMl 
measures a lifetime Tp associated with the decay of pho- 
toexcited electron population. Systematically compar¬ 
ing Ts{e) and rp(e) may provide new insights into the 
underlying scattering mechanisms. To understand the 
relation between these two lifetimes, one needs to ob¬ 
tain the energy-resolved population lifetime Tp(e) and di¬ 
rectly compare with rs(e). Femtosecond time-resolved 
ARPES (trARPES) provides this capability Sev¬ 

eral trARPES studies have investigated the relaxation of 
photo-excited electrons in cuprates [IMi. Yet, so far, 
no energy-resolved lifetimes have been extracted from the 
population dynamics in cuprates. 

In this Letter, we employ trARPES and ARPES 


with high energy resolution to perform a detailed com¬ 
parison between rs(e) and Tp(e) in optimally doped 
Bi2Sr2Cao.92Yo.o8Cu208+5 (OP Bi-2212, T, = 96 K) 
along the nodal direction [25|- At 20 K, rp(e) extracted 
from trARPES decreases with increasing excitation den¬ 
sities below a characteristic energy of ~ 60 meV, yet the 
trend is reversed above this energy. At first glance, this 
characteristic energy seems to agree with the mode ener¬ 
gies as identified by ARPES measurements of Ts{e), but 
the absolute values for rs(e) and Tp(e) are different by 
1 2 orders of magnitude. This di sparity also existed 


in studies on graphite and graphene |26l - l28j | . We demon¬ 


strate that Ts(e) and rp(e) reflect different aspects of elec¬ 
tron scattering phenomena and that processes beyond 
electron-phonon interactions contribute to the disparity. 
The understanding of this disparity is of importance to 
future trARPES experiments on all materials. 


Our trARPES setup is based on a Ti-Sapphire regener¬ 
ative amplifier operating at a repetition rate of 800 kHz. 
1.5 eV infrared pump pulses excite the sample; 6 eV ul¬ 
traviolet probe pulses generate photoelectrons which are 
collected by a Scienta R4000 analyzer. High quality sin¬ 
gle crystals of OP Bi-2212 are cleaved in ultrahigh 
vacuum with a pressure < 7 x 10“^^ Torr. Typical en¬ 
ergy, momentum, and time resolutions for the trARPES 
setup are 22 meV, 0.001 A“^, and 100 fs, respectively. 
Our ARPES measurement is performed at Beamline 5- 
4 of Stanford Synchrotron Radiation Lightsource. 7 eV 
synchrotron light generates photoelectrons which are col¬ 
lected by a Scienta R4000 analyzer. The sample prepa- 
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FIG. 1. Energy-resolved population decay analysis, (a) Scheme of a pump-probe photoemission experiment, (b) Nodal cut 
obtained by 6 eV photoemission at 20 K. Boxes mark the windows for momentum integration and 8 meV energy binning used 
for the population decay analysis. (c)~(e) Population dynamics normalized by peak intensities at 20 K for incident excitation 
densities (fiuences) of 10, 27, and 60 ytJ.cm~^, respectively. (g)~(i) Population dynamics at 120 K for the same set of excitation 
densities, (f, j) Population decay rates extracted by fitting the transients with a Gaussian-convolved exponential decay at initial 
delays (A/(t) > 20%A/niax). The fitting curves (thick black lines) are overlaid on the population dynamics. 


ration procedure is the same as that for the trARPES 
measurement. The combined energy resolution is 3 meV. 

We study the cut along the nodal direction in Bi-2212 
as shown in Fig.[ljb). A clear kink in the band dispersion 
^ 70 meV below the Fermi level E-p can be identihed. 
This indicates a strong bosonic cormling and agrees with 
previous ARPES investigations [J-lSl- 

trARPES measurements are also performed on the 
nodal cut. For the following analysis and discussion, we 
focus on the photoexcited electron population above Ep 
[^ . We integrate over the whole momentum range of 
the cut as illustrated in Fig. [TKb) and subtract the sig¬ 
nal before pumping. Energy-resolved transient electron 
populations are obtained by binning the data into 8 meV 
energy intervals and plotting as a function of pump-probe 
delay. The population dynamics normalized by their 
peak intensities are shown in Fig. [IKc)-(e) for incident 
fiuences of 10, 27, and 60 /iJ.cm”^ at T = 20 K. The 
population dynamics for the same set of fiuences at 120 K 
are displayed in Fig. [T](g)-(i). These population dynam¬ 
ics are fitted with an exponential decay convolved with a 
Gaussian function for initial delays as defined by a 20% 
intensity cutoff [IlElii. This yields the population 
decay rates as a function of energy, which are displayed 


in Fig. [TJf) and (j). The fit results are insensitive to the 
choice of cutoff values [^ . 

We first examine the energy dependence of the pop¬ 
ulation dynamics. A pronounced energy dependence is 
most clearly observed for the 10 /xJ.cm”^ data at 20 K 
(Fig.IHc)). Notably, an abrupt change of the population 
dynamics occurs near 60 meV. Below this energy the 
populations increase for ~ 0.5 ps before reaching their 
maxima, and live as long as a few ps; above this en¬ 
ergy the populations reach their maxima near time zero, 
and decay within a few 100 fs. Consistently, the corre¬ 
sponding decay rate in Fig. [IJf) displays a pronounced 
increase by one order of magnitude near 60 meV. The 
energy dependence for the same fluence at 120 K is less 
drastic (Fig.[T](g)). Populations at all energies reach their 
maxima near time zero, and decay in a few 100 fs. Never¬ 
theless, an abrupt increase near 60 ~ 80 meV is observed 
in the decay rate (Fig. [T](j)). 

We then study the excitation-density dependence of 
the population dynamics. The excitation density is char¬ 
acterized by the pump fluence which specifies the incident 
energy per unit area. At 20 K, the rising edges become 
gradually less delayed as the fluence increases. More¬ 
over, the extracted decay rates display a pivoting behav- 
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FIG. 2. Comparison between ImE (black lines) obtained by 
ARPES below Ef and population decay rates (solid circles) 
obtained by trARPES above Ef- The method of extracting 
ImE from the ARPES data is described in the Supplemen¬ 
tal Material [^ . Note that the population decay rates are 
magnified by a factor of 10 for visualization purposes. 


ior when tuning the pump fluence (Fig. [IJf) and [21(a)). 
Below 60 meV, the decay rate increases with increasing 
fluence. Above 60 meV, the decay rate weakly decreases. 
This pivoting is much weaker at 120 K, where the decay 
rates are approximately fluence independent (Fig. [ijj)). 

We summarize the three key observations: (i) a ^ 
0.5 ps rising edge at the lowest fluence and temperature; 
(ii) an abrupt increase in decay rates near 60 ^ 80 meV 
at the lowest fluence; (iii) a pivoting behavior for the 
rate curves when tuning the pump fluence. This com¬ 
plex energy and fluence dependences establish a multi¬ 
dimensional constraint for a microscopic understanding. 

A delayed rising edge is usually attributed to cascade 
processes which fill low-energy states using high-energy 
electrons [2^,[3lj. For cuprates, the d-wave gap indicates 
that the available low-energy electronic states are only 
near the node, which provides an additional constraint 
for electron accumulation. This constraint should be 
lifted if the superconducting gap is melted either by rais¬ 
ing the equilibrium temperature above Tc, or by pump¬ 
ing beyond the fluence of ISuJ.cm”^ which is required 
to transiently melt the gap [l9!,l20j. Indeed, Fig. |T] shows 
that the rising edges of the population dynamics are sig¬ 
nificantly shortened in both situations. 

The abrupt change and the pivoting point of Tp occur 
at the same energy ~ 60 meV, which is reminiscent of 
mode energies revealed by ARPES measurements of Ts- 
The pivoting behavior of Tp due to electron-phonon cou¬ 
pling has been theoretically predicted [s^. This model 
further predicts that and Tp converge in the zero exci¬ 
tation limit, yet we find a significant quantitative differ¬ 
ence (Fig. 12). In particular, if we compare h/{2Tp) and 
ImE below 50 meV, instead of converging fi/{2Tp) devi¬ 
ates further from ImE when lowering the fluence. While 
ImE extracted from ARPES depends on photon energy, 
it is always on the same order of magnitude 0|. In 


contrast, the discrepancy between Tp from tr ARPES and 
Ts from ARPES is 1^2 orders of magnitude. 

While one may relate this discrepancy to unique scat¬ 
tering properties in cuprates, we notice that it is indepen¬ 
dent of particular material systems. As shown in Table U 
the discrepancy between Tij (2Tn) and ImE exists also in 
graphite and graphene |26l428l| . 


TABLE I. Comparison of scattering rates at the apparent 
mode energies (Ho) obtained by ARPES and trARPES on 
several materials 


Material 

Ho (meV) 

ImE (meV) 

V(2- 

rp) (meV) 

OP Bi2212 

70 

33 


1.8 

Graphite [26, 27} 

200 

170 


1.6 

Graphene [2^ 

200 

118 


0.37 


The generality of this discrepancy poses a challenge for 
connecting time-resolved experiments to single-particle 
scattering properties. In the following we discuss the 
conceptual difference between Tp and Tg, and survey rep¬ 
resentative scattering channels of importance to all ma¬ 
terials. 

Collision integrals [H, 33-3^ provide a general formal¬ 
ism for electron scattering processes. 


/ dk' 

^IFkqk/(ek0[l-/(ek)] (1) 

Ek denotes the electronic state at momentum k on a 
band dispersion, /(ck) is the corresponding occupation. 
Wk.k' stands for the probability of scattering from k to 
k'. The two integrals in Eqn. [T] represent the emptying 
processes from the state at Ek to other states, and the 
filling processes from other states back to the state at Ek. 

In equilibrium, all the scattering processes reach a 
detailed balance such that df{€-k)/dt = 0. In non¬ 
equilibrium and for perturbative excitations, /(ck) and 
/(Ek') are transiently changed to /(ck) + Sf{ek) and 
/(ek') + ^/(ek')- Eqn.|T]is expanded to first order in S f to 
describe the evolution of photoexcited electrons [s^, . 

According to the definition of Ts(Ek), it is associated with 
excitations that only change /(ck), and hence Sf{e\^') = 0 
in the Taylor expansion. This occurs in an ARPES mea¬ 
surement where the incident photon probes the same 
photo-hole as it excites 0. Using this concept, we derive 
a general expression for Ts(Ek). 

1 _ 1 di5/(Ek) 

'rs(ek) <5/(ek) dt 
f dk' 

= J ^{VFk.k'[l-/(ek')] + IUk',k/(ek')} 

( 2 ) 
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FIG. 3. Illustration of single-particle scattering and popula¬ 
tion evolution. (a)'~(c) Single-particle scattering via electron- 
impurity, electron-electron, and electron-phonon interactions. 
The scattering rate at arises from the non-equilibrium dis¬ 
tribution at ek. (d)'~(f) Population evolution via electron- 
impurity, electron-electron, and electron-phonon interactions. 
The population decay rate at ck arises from non-equilibrium 
distributions at tk and other states at Ck'. Red boxes denote 
the integration window as used in our data analysis. 


On the other hand, in time-resolved measurements, the 
pump pulse can excite electrons into many different states 
independent of the probe pulse. Both (5/(ek) and (5/(ek') 
are in general non-zero, which leads to Eqn. [3]for Tp(ek). 


dk‘ 


Tp(ek) 'rs(ek) 


(27r) 


3 {Wk,k'/(ek) 


+l^k',k[l-/(ek)]} 


<5/(ek) 


(3) 


Eqn. [3] demonstrates the general distinction between 
Tg and Tp without specifying the scattering mechanism. 
This result is consistent with a number of theoreti¬ 
cal investigations which numerically or analytically 
[3 III solve the collision integrals. We emphasize that 
when (5/(ek') and 5/(ek) individually approach zero in 
the zero excitation limit, their ratio can be nonzero. 
Therefore, the distinction between Tg and Tp is funda¬ 
mental. In the following we illustrate this distinction by 
surveying a few common scattering channels. 

Electron-impurity scattering. This channel is often 
considered to be elastic and energy-independent [Ij. 
Since Wk.k' = lEkqk, Eqn. [5] and [3] are significantly 
simplified: 


' s,e-imp 


(^k) 


dk' 

(27r)3 


Wk,k' 


(4) 


.(ek) 


/ 


dk' 

(27r)' 


dEk,k'[l- 


^/(£k0 . 

^/(ek) ^ 


(5) 


Notably, the electron-impurity scattering makes a 
nonzero contribution to ImE (Eig. [2Ia)), but less to 
7i/(2Tp). The latter contribution completely vanishes for 
isotropic excitation ((5/(ek')/'^/(£k) = 1), which results 
from symmetric scattering between k and k' (Fig. [3l^d)). 
This is consistent with conclusions reached by spherical 
harmonic decomposition of /(ck) [3 38|. Removing a 


constant from the experimental ImS may let us reconcile 
the discrepancy between ImE and 7i/(2rp) for low ener¬ 
gies (Fig. [2]). However, the relative increase of ImE(e) 
from e = 14 meV to e = 70 meV is about 20 meV, which 
is still ^ 10 times larger than the counterpart for fij (2Tp). 

Electron-electron scattering. A photoexcited electron 
can interact with another electron to redistribute energy 
and momentum (Fig. [3Kb)). Time-resolved photoemis¬ 
sion experiments on metals have demonstrated that the 
population lifetimes for electrons > 0.5 eV above Ey are 
dominated by electron-electron scattering, and can be 


related to the single-particle lifetimes [39|, . However 


for low-energy electrons the filling processes from below 
Ey become important (Fig. |3Ke)), which slow down the 
population decay. 

Electron-phonon scattering. A photoexcited electron 
can dissipate its energy by interacting with phonons of 
energy H (Fig. EKc)). ARPES experiments [3^ have 
concluded that this is an important scattering channel 
in cuprates for low-energy electrons < 0.1 eV. The dis¬ 
tinction between Tp and Tg in Eqn. |3| still applies. How¬ 
ever, electron-phonon scattering is a special case where 
this distinction can be mitigated in the low-excitation, 
low-temperature limit 3 32, ^|. Scattering processes 
associated with the photoexcited populations i5/(ek') at 
Ck' = Ck + and Ck' = Ck — affect the population 
at Ck in the form of phonon emission and absorption, 
respectively. In the low-temperature limit, the phonon 
absorption processes are negligible due to the vanishing 
Bose-Einstein distribution. In the low-fluence limit, the 
phonon emission processes from Ck' = Ck + to Ck are 
also negligible since the population at Ck + H can be ne¬ 
glected compared to that at e^. 

The three cases above demonstrate that Tp is in general 
distinct from Tg. Electron-phonon scattering is a special 
case where Tp converges to Tg in the low-temperature, 
low-fluence limit. Notably, when adding more scattering 
channels the total single-particle scattering rate increases 
additively. Yet this does not apply to the population de¬ 
cay rate. The main difference, as shown in Eqn.jSl is that 
certain scattering channels can contribute negatively to 
the total population decay rate due to filling processes, 
and slow down the overall population decay. Therefore, 
although we have a qualitative understanding for each 
individual scattering channel, understanding their inter- 
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play in the population decay is nontrivial. For example, 
while electron-phonon coupling alone can give rise to a 
pivoting point of Tp, the origin of the pivoting point in 
our data remains unclear. 

The discrepancy between Tg and Tp in our data sug¬ 
gests that scattering channels beyond electron-phonon 
scattering have an appreciable contribution. Electron- 
impurity and electron-electron scatterings are likely to 
contribute. In the superconducting state Wk,k' includes 
pair breaking and reformiM processes, and a supercon¬ 
ducting gap opens in Ck The dynamics in the su¬ 

perconducting state can be addressed by the phenomeno¬ 
logical Rothwarf-Taylor model, which describes electron 
pai ring and boson-induced pair breaking processes [13- 
Q, lldl . [13 . [l3. [ 43 . Hi ■ A number of trR 12 - 13| and 
trARPES [is], [l3| experiments have demonstrated that 
the energy-integrated population lifetimes in cuprates de¬ 
pend on the pump fluence, which is also the case in our 
data [13]. This has been interpreted as dynamics not 
governed by the boson bottleneck 12 However, due 
to the nonlinear nature of the Rothwarf-Taylor model, 
the fluence-dependent decay rates may be present in the 
bottleneck limit as well [la, [44 1. To make direct com¬ 
parison with our energy-resolved rates, a more complete 
model incorporating the energy-dependence of the pair¬ 
ing processes will be required. 

The unique energy- and fluence-dependences of our 
data serve as a basis for understanding scattering phe¬ 
nomena in cuprates. The delayed rising edges at low 
temperature reflect the anisotropic gap. The discrepancy 
between and Tp suggests scattering channels beyond 
electron-phonon interactions. Within the formalism we 
utilize, our data provides a strong constraint for future 
theories exploring the specific form of IFk.k'- Our work 
thus points to a new route to unveil the dominant scatter¬ 
ing channels behind high-temperature superconductivity. 
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